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Abstract 

We study a three-dimensional system of particles interacting via spherically- 
symmetric pair potentials consisting of several discontinuous steps. We show that at 
certain values of the parameters desribing the potential, the system has three first- 
order phase transitions between fluids of different densities ending in three critical 
points. 



Much attention has been focused on the topic of liquid-liquid phase transi- 
tions [1]. Recently a liquid-liquid phase transition was found experimentally 
in phosphorus [2]. Liquid-liquid phase transitions also can exist in other ma- 
terials with tetrahedral symmetry such as carbon [3,4], silica, [5] and silicon 
[6,7]. 

The possibility of a sharp amorphous-amorphous transition between low den- 
sity amorphous ice (LDA) and high density amorphous ice (HDA) has been 
known for a long time [8]. Computer simulations of the ST2 water model 
[9] suggest that the transition line between LDA and HDA can be extended 
to higher temperatures above the glass transition, where it would become 
a first order liquid-liquid phase transition line ending at a critical point at 
about 200K [10]. However, in this temperature range, water immediately crys- 
tallizes and this phase transition cannot be observed in direct experiments. 
More recently, a distinct amorphous phase of water, called very high-density 
amorphous ice, has been observed experimentally [11-13]. New simulations of 
the ST2 model [14] in the restricted NPT ensemble suggest that instead of 
one previously observed liquid-liquid phase transition, there might be three 
liquid-liquid phase transition lines between four liquids with different den- 
sities roughly corresponding to the densities of the three known amorphous 
ices. It has been suggested [15] that systems with rich polymorphic solid phase 
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diagrams may have several liquid phases with local particle arrangements cor- 
responding to various crystalline structures. 

We propose here a simple model with at least two liquid-liquid phase transi- 
tions in addition to the known liquid-gas phase transition. This model is based 
on the double step soft-core potential model previously used to study systems 
with two fluid- fluid critical points [16-19]. Two fluid-fluid critical points were 
also observed by Monte Carlo simulations of a similar model with a continuous 
soft-core potential [20]. 

We perform constant volume V collision-driven molecular dynamic simulations 
[17,18,21] on a system of N = 1728 spherically symmetric particles interacting 
via the discontinuous pair potential 
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where U(r) is a potential energy of a pair of particles at distance r, Ua > 
is the energy of attraction, U\ = 3Ua, U 2 = Ua, &i = 1.4a, b 2 = 2.1a, and 
c = 2.8a. (see Figure 1, inset). We keep temperature T constant using a mod- 
ified Berendsen method and compute pressure P the same way as in [16,17], 
averaging it over 2 • 10 4 (ma 2 /Ua) 1 / 2 time units, where m is particle mass. 
Throughout this paper, we measure density, pressure, and temperature in 
units of a -3 , f/^a -3 , and U a/&b, respectively, where ks is the Boltzmann con- 
stant. We investigate several state points for systems with different densities 
p = N/V along several isotherms (see Fig. 1) and find three distinct sets of 
van der Waals loops, each ending in a specific critical point characterized by 
an isotherm with an inflection point. 

Note that the isotherm for T = 0.60 almost entirely corresponds to a stretched 
metastable liquid with P < 0. In addition, at this temperature, the system 
spontaneously crystallizes for p > 0.23. Nevertheless, we are able to find equi- 
librium liquid pressures by averaging pressure before the crystallization onset, 
using the same technique as in [17]. We quench five independent configurations 
equilibrated at T=0.9 to T=0.6, skip initial 800 time units (the correlation 
time at this temperature), and average the pressure until the onset of crys- 
tallization, indicated by the sharp decrease of potential energy. The average 
existence time (t) of the metastable liquid is (t) = 1.5 • 10 4 time units at 
p = 0.23 , (t) = 4.5 • 10 3 at p = 0.25, (t) = 1.5 • 10 3 at p = 0.27, and (t) = 10 3 
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Fig. 1. Isotherms and critical points for the system of spherically symmetric particles 
interacting via the potential shown as an inset. 

at p = 0.30. The error bars of the pressure are assumed to be equal to the 
standard deviation of the pressure measurements for five independent runs. 

Accordingly, we find a liquid-gas critical point C 1 at Ti = 1.10, Pi = 0.014, 
Pi = 0.054, a stable high density critical point C 2 at T 2 = 0.75, P 2 = 0.0056, 
P2 = 0.152, and a third very high density critical point C3 at T 3 = 0.60, 
P3 = —0.014, P3 = 0.245. The third critical point is metastable with respect 
to gas-crystal coexistence, and by varying the parameters of the potential [19], 
one should be able to shift the position of the third critical point to a stable 
region of the phase diagram. 

Figure 2 shows radial distribution functions for the three liquids at temper- 
ature T = 0.59, corresponding to the low density liquid (LDL) with density 
0.098a -3 , high density liquid (HDL) with density 0.186a~ 3 , very high density 
liquid (VHDL) with density 0.272a~ 3 and a crystal with density 0.296a -3 . One 
can see a dramatic difference between the structures of LDL and HDL. 

Integrating radial distribution function for LDL, we find that 11 neighbors 
stay in the attractive well, almost not penetrating into soft cores. 

In contrast, in HDL, particles penetrate into a larger soft core. The average 
number of neighbors in the smaller soft core is 0.6. The average number of 
neighbors in the larger soft core is 2.7. The number of particles in the attractive 
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Fig. 2. Radial distribution functions for a system shown in Fig. 1 at temperature 
T = 0.59 at four different densities, corresponding to three different liquids and a 
crystal. 

well also increases to 17. 

The difference between the structures of HDL and VHDL is less pronounced. 
Nevertheless, one can see the development of an extra peak for VHDL at a 
distance 3.5a. In addition, the number of particles in the smaller and larger 
soft cores increases to 1 and 5, respectively, and the number of particles in the 
attractive well increases to 23. 

The choice of parameters is motivated by our preliminary studies of sim- 
pler systems with U\ — oo, which have a gas-liquid critical point C\ ( T — 
1.06, P = 0.013, p = 0.052) and a stable liquid-liquid critical point C 2 
(T = 0.67, P = 0.059, p = 0.144 ) (see Fig. 3). The potential of this system 
is the soft-core potential [16], with b\ playing the role of the hard core and 
b 2 playing the role of the soft core. Note that the positions of these critical 
points are very similar to those of C\ and C 2 in the three-step potential. Thus, 
adding a second soft core a < b x < b 2 to the potential, creates a third critical 
point at higher densities without a significant change in the position of the 
two original critical points. 

One can speculate that adding more steps to a potential can create more 
critical points at different densities corresponding to the interparticle distances 
of these steps if the parameters of these steps are carefully selected. 
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Fig. 3. Isotherms and critical points for temperatures T =1.5, 1.45, 1.4, 1.35, 1.30, 
1.25, 1.20, 1.15, 1.10, 1.05, 1.00, 0.95, 0.90, 0.85, 0.80, 0.75, 0.70, 0.65, 0.60 for a sim- 
plified potential shown in the inset with one soft core hi and a hard core b\ = 1.4a. 
Isotherms for T > 1.00 are computed only for small densities, while isotherms for 
T < 1.00 are computed only for large densities. The isotherm corresponding to 
T = 1.00 is computed for the entire range of densities. The systems for p > 0.23 
spontaneously crystallize at T < 0.75 into crystalline structures similar to to the 
crystalline structure of the system with two soft cores shown in Fig. 2. 



In summary, we show that systems interacting via spherically symmetric po- 
tentials with three characteristic repulsive distances playing the role of soft 
cores may have three liquid phases of three increasing densities characterized 
by the penetration of particles into a soft core of smaller diameter. 
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